AREAS OF CONTACT AND PRESSURE DISTRIBUTION 


IN BOLTED JOINTS 


H.H. Gould 
B . B . Mikic 


Final Technical Report 
Prepared for 


George C. Marshall Space Flight Center 
Marshall Space Flight Center 
Alabama 35812 


under 


Contract NAS 8-24867 


June 1970 


DSR Project 71821-68 


Engineering Projects Laboratory 
Department of Mechanical Engineering 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 




- 2 - 


ABSTRACT 


When two plates are bolted (or riveted) together these will be in 
contact in the immediate vicinity of the bolt heads and separated 
beyond it. The pressure distribution and size of the contact zone is 
of considerable interest in the study of heat transfer across bolted 
joints . 


The pressure distributions in the contact zones and the radii at 
which flat and smooth axisymmetric , linear elastic plates will separate 
were computed for several thicknesses as a function of the configuration 
of the bolt load by the finite element method. The radii of separation 
were also measured by two experimental methods. One method employed 
autoradiographic techniques. The other method measured the polished 
area around the bolt hole of the plates caused- by sliding under load 
in the contact zone. The sliding was produced by rotating one plate 
of a mated pair relative to the other plate with the bolt force acting. 


The computational and experimental results are in agreement and 
these yield smaller zones of contact than indicated by the literature. 
It is shown that the discrepancy is due to an assumption made in the 
previous analyses. 

In addition to the above results this report contains the finite 
element and heat transfer computer programs used in this study. 
Instructions for the use of these programs are also included. 
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NOMENCLATURE 


A, 

D 

E 

G 

h 

c 

H 

k, 

P, 


B, C 


h. 



r 


R 

o 

u,w 

X 


X 

c 

y 

y' 

z 


radii 

thickness 

modulus of elasticity 
shear modulus 

heat transfer coefficients 
hardness 

thermal conductivities 

pressure 

coordinate 

radius of separation 

displacement in r and z 

coordinate 

length of contact 

coordinate 

slope 

coordinate 



deflection 

dilation 

strains 

standard deviations 


directions 


stresses 
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X , \i 

Lame’s constants 

V 

Poisson’s ratio 

T 

shear stress 

0 

angle 

Subscripts 

r 

radial direction 

t 

tangential direction 

z 

z-direction 
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chapter I 
INTRODUCTION 

When two plates are bolted (or riveted) together, these will be 
in contact in the immediate vicinity of the bolt heads and separated 
beyond it. The pressure distribution in the contact area and the 
separation of the plates is of considerable interest in the study of 
heat transfer across joints. Cooper, Mikic and Yovanovich [1] show 
that with assumed Gaussian distribution of surface heights, the micro- 
scopic contact conductance is related to the interface pressure, 
surface characteristics and the hardness of the softer material in 


h 

c 


1.45 


tan 6 

a 


p 0 . 98 5 


( 1 . 1 ) 


where 


k 


2k., k^ 


( 1 . 2 ) 


and k-, and k„ represent the thermal conductivities of two bodies in 
contact; o is the combined standard deviation for the two surfaces 
which can be expressed as 


/ 2 , 2 , 
= (a^ + O^) 


O 


1/2 


(1.3) 
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where and are the individual standard deviation of height for 

the respective surfaces; tan 0 is the mean of the absolute value of 
slope for the combined profile and it is related, for normal distribu- 
tion of slope, to the individual mean of absolute values of slopes as 


tan 0 = 


2 2 
(tan 6^ + tan G^) 


1/2 


(1.4) 


where 


tan 0 . 
1 


lim 

L 


OO 


1 

L 


.L 

|y^| dx; i = 1,2 (1.5) 


and y' is the slope of the respective surface profiles; P represents 
the local interface pressure; and H is the hardness of the softer 
material. 

Relation (1.1), as written above, is applicable for contact in a 
vacuum. One can modify the expression by simply adding to it 


_ conductivity of interstitial fluid 
f ~ average distance between the surfaces 


( 1 . 6 ) 


in order to account approximately for the presence of the interstitial 
fluid. 

All parameters in relation (1.1), except for the pressure, are 
functions of the material and geometry and can be easily obtained. The 
determination of the pressure distribution and the extent of the contact 
area between two plates present both mathematical and experimental 
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difficulties. From the mathematical point of view, the difficulty stems 
from the fact that the theory of elasticity will yield a three dimensional 
(axisymmetric) problem with mixed boundary conditions. Experimentally, 
the discrimination between contact and gaps of the order of millionths 
of an inch is required. 

Roetscher [2] proposed in 1927, a rule of thumb that the pressure 
distribution of two bolted plates, Fig.l, is limited to the two frustums 
of the cones with a half cone angle of 45 degrees as shown in Fig. 2 and 
that at any level within the cone the pressure is constant. Also, for 
symmetric plates, according to Roetscher, separation will occur at the 
circle which is defined by the contact plane and the 45 degree truncated 
cone emanating from the outer radius of the bolt head. 

Since 1961 Femlund [3], Greenwood [4] and Lardner [5] among others 
reported solutions based upon the theory of elasticity. Although their 
solutions also yield separation radii at approximately 45 degrees as in 
Roetscher’ s rule, their solutions yield a much more reasonable pressure 
distribution as compared to Roetscher’ s constant pressure at each level 
of the frustrum. These investigators have made use of the Hankel 
transform method demonstrated by Sneddon [6] in his solution for the 
elastic stresses produced in a thick plate of infinite radius by the 
application of pressure to its free surfaces. The basic assumption in 
their approach is that two bolted plates can be represented by a single 
plate of the same thickness as the combined thickness of the two plates 
under the same external loading. It then follows tnat the z-stress 
distribution at the parting plane can be approximated by the z-stress 
distribution in the same plane of the single plate. It also follows 
that separation will occur at the smallest radius in that plane for which 
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the z-stress is tensile. In the case of two plates of equal thickness 
the stress at the midplane of the equivalent single plate is the 

stress of interest. 

Feinlund [3], for example, used the method of superposition in the 
sequence shown in Figs. 3(a) to 3(c) to obtain annular loading. Then 
by superposition of shear and radial stresses at radius A, Figs. 3(d) 
and 3 (e) , opposite in sign of those due to the annular loading at the 
free surfaces, Fernlund obtained the solution for a single plate with a 
hole under annular loading (Fig. 3(f)). 

Experimental work in this area included Bradley’s [7] measurements 
of the stress field by three dimensional photoelasticity techniques, and 
the use of introducing pressurized oil at various radii in the contact 
zone and measuring the pressure at which oil leaks out from the joint 
[3,8]. Both of these experimental methods have uncertainties as indi- 
cated by the authors. 

Because of the cumbersomness of the Hankel transform solution and 
experimental difficulties, the body of work in this area has been very 
limited and definite verification of analytical results by experiment is 
not cited in the literature. 

The research described in the succeeding chapters was undertaken 
with the following primary objectives: 

a) To provide a method of solution for the case of two bolted 
plates without the simplifying assumption of the single 
plate substitution. 

b) To devise a test to validate the two plate analysis. 

c) To test the validity of the single plate substitution. 
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A finite element computer program has been assembled for the 
analytical solution of two-plate problems. Experiments have been 
performed to verify the analytical results. Since in heat transfer 
calculations the extent of the radius of contact is of primary 
importance, and since by restricting the experimental effort to the 
verification of only this parameter, (rather than the verification of 
the entire pressure distribution,) many experimental uncertainties 
should be eliminated, the experiments were designed only for the deter- 
mination of the contact area. 

Agreement between analysis and experiment was obtained and the 
results show that the single plate substitution is not justified and 
the 45 degree rule is not valid for the flat and smooth surfaces 
studied. 
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Chapter II 
ANALYSIS 

A. Problem Statement 

The objective of the anlaysis was to solve the linear elasticity 
problem of two plates in contact defined mathematically by the following 
equations for each plate: 

The equations of equilibrium 


/ (ra ) - 
or r 

^ , 9t 

'’t + ^ 

= 0 

(2.1) 


^ (r a ) 
9z z 

= 0 


where i = T 

rz zr 

= T 

and T = T = T 

rt tr zt 

= T =0. 

tz 

The stress - strain 

relations, 

using standard notation 

for stress and 

strain, 




a = Ae 
r 

+ 2 ye 

r 




+ 2 ye^ 


(2.2) 

II 

00 

+ 2 ye 

z 



T = 2 ye 

rz 
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where X and y are Lame’s constants and 


2 GV 

1 - 2v 


y = G 


(2.3) 


if G is the modulus of elasticity in shear and V is Poisson's ratio 
and e the volume expansion is defined by 


e = 


3u u 9w 

^ F ^ 


(2.4) 


where u is the displacement in the radial direction and w is the 
displacement in the axial direction. 

The strain - displacement relations 


e 

r 


9u 

9r 


e 


t 


u 

r 


S 9z 

1.9u , 9w V 
^rz “ 2^3z 3r ^ 


(2.5) 


The above equations can be combined to yield the equilibrium 
equations in terms of displacements 


V^u 


u , 1 9 e 

2 1 - 2v 9r 

r 


0 


V^w + 


1 - 2v 


9s 

9z 


( 2 . 6 ) 


0 
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The applicable boundary conditions are (see Fig. 11) 


(A,z) = 

r r 

(A,z) = 

o^^hc.z) = (C,z) = 

r r 

T^^hc.z) = (C.z) = 

(r, -D^) 

T^^^hr.O) = (r.O), 

ofhr.Di) = (r, -D^) 

T*^^\r,0) = x'^^hr.O) = 

a^^hr.O) = O^^hr.O), 
a^^hr.O) = af^hr.O) = 

z z 

afhr.D^) = afhr, -D^) 

w^^\r,0) = w^^^(r,0). 


0 

0 

0 

0 

= 0 , 

A < r < R 

_ __ o 

= 0, b j< r £ C 

(2.7) 

0, R < r < C 

o— — 

A < r < R 
— — o 

0, R < r < R 

o— — 

= P(r) , A r B 

A < r < R 



A 


A 
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Inspection of the above equations shows that the above constitutes 
a mixed boundary value problem and the most appropriate technique for 
solution is the finite element method. 

B, Method of Analysis 

A finite element computer program was assembled for the analytical 
solution of bolted plates. Descriptions of the finite element method 
are given in references [9,10], but for completeness, an outline of the 
mathematical formulation for this case is presented in Appendix A. A 
listing of the cumputer program and instructions for its use may be 
found in Appendix B. Appendix C contains user’s instructions and a 
listing of the finite element program modified to include thermal strains. 

As in the previous work axial S 3 nnmetry and isotropic linear elastic 
material behavior were assumed. However, the computer programs accom- 
modate plates with different material properties in a bolted pair. 

The basic concept of the finite element method is that a body may 
be considered to be an assemblage of individual elements. The body then 
consists of a finite number of such elements interconnected at a finite 
number of nodal points or nodal circles. The finite character of the 
structural connectivity makes it possible to obtain a solution by means 
of simultaneous algebraic equations. When the problem, as is the case 
here, is expressed in a cylindrical coordinate system and in the 
presence of axial symmetry in geometry and load, tangential displacements 
do not exist, and the three-dimensional annular ring finite element is 
then reduced to the characteristics of a two-dimensional finite element. 
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The analysis consists of (a) structural idealization, (b) evaluation 
of the element properties, and (c) structural analysis of the assem- 
blage of the elements. Items (b) and (c) are covered in the appendices 
and in the references quoted. The structural idealization and the 
criteria for acceptable solutions will be described in this chapter. 

Fig. 4(a) shows two circular plates in contact under arbitrary 
axisymmetric loading. The plates are subdivided into a number of 
annular ring elements which are defined by the corner nodal circles 
(or node points when represented in a plane) as shown in Figs. 4(b) and 
4(c). Unlike the cases described in Chapter I, which have been solved 
by the Hankel transform method, all plates solved by the finite element 
method have finite radii. The cross sections, of each annular ring 
element is either a general quadrilateral or triangle. To improve 
accuracy smaller elements are used in zones where rapid variations in 
stress are anticipated than in zones of constant stress; thus the 
different size elements shown in Fig. 4(b). (However, the total number 
of elements allowable are subject to computer capacity.) 

Figure 4(b) shows the two plates in contact for the radial distance 

X and separated beyond it. It is to be noted that the nodal points 
c 

on the parting line and within the length of contact are common to 

elements in both plates. The other elements adjacent to the parting 
line on each plate are separated from their corresponding elements in 
the mating plate and these elements have no common nodal points. 
Physically, it is equivalent to the welding together of the two plates 
in the contact zone. Mathematically, we are imposing the condition that 
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in the contact zone the displacements in the z and r directions be 
identical for both plates. In the case of bolted plates of equal 
thickness, i.e. in the presence of symmetry about the parting plane, 
these conditions apply exactly. Furthermore, because of this symmetry, 
one needs to analyze only one plate, as shown in Fig. 5(b), with the 
imposed boundary conditions on the contact zone of zero displacement in 
the z-direction and freedom to displace in the r-direction. It can 
also be observed that the solution of two plates with symmetry about the 
parting plane is equivalent to the solution of one of these plates under 
the same loading conditions, but resting on a frictionless infinitely 
rigid plane. Also, under the above conditions the shear stress in the 
contact zone is identically zero. 

In the case of bolted plates of unequal thickness the model includes 
both plates as shown in Fig. 5(c). This model is an approximation 
because, in general, two plates of mequal thickness do not have the same 
displacement in the r-direction on the contact surface. The solution 
yields, therefore, a shearing stress distribution in the contact zone. 

The solution, however, should be exactly compatible with the physical 
model if the frictional forces in the joint prevent sliding. 

The critical aspect of the approach used herein is the determination 
of the largest nodal circle on the parting plane which is common to an 
element on each plate. This nodal circle defines the contact zone and 
the radius, R^, at which separation occurs. 

The output of the finite element computer program includes the 
displacement of each node in the r and z directions and the average 
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a , a , and T stresses for each element, 

z r t rz 

The computation is iterative and the objective is to achieve the 

lowest possible compressive stress in the outermost elements 

bordering the contact zone. Unacceptable solutions are shown in Fig. 6(a) 

and 6(b). If R for a given external load distribution is too small, 
o 

then the solution will show that the two plates intersect (Fig. 6(a)). 

On the other hand, if is assumed too large, the solution will show 

that the outer portion of the contact zone sustains a tensile cr stress 

z 

(Fig. 6(b)). Neither of these two situations is physically feasible. 

In general, the procedure employed was to commence the iterations with 

a value for R which would yield a tensile a stress in the outer 
o z 

elements adjacent to the contact zone and then move R inward. The 
iteration ended as soon as no tensile stress was present at the 

contact zone. For example, for the case shown in Fig. 5(b), if the 
stress for the element in the last row and to the left of the last roller 
is tensile, then the following iteration will proceed without the last 
roller. Thus, the resolution is one nodal interval. Finer resolution 
can be obtained by reducing the interval between nodal circles by 
introducing more elements or shifting the grid locally. The same 
criteria apply to the model shown in Fig. 5(c). 

In the finite element analysis of the Fernlund (3) model, i.e. 
single plate with external loads at the faces z = + D no iteration is 
required and the rollers shown in Fig. 5(c) would extend to the outer 
radius of the plate. (Although Fernlund’ s computations are based on 
infinite plates, computations show that there is no distinction between 
infinite plates and plates of radius greater than five times of the outer 
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radius, B, of the load. See Fig. 5(a). 

Convergence was tested by subdividing elements further, with nodal 
points in the coarser grid remaining nodal points in the finer grid. 
Changing the mesh from 180 elements to 360 elements have shown no 
improvement in accuracy. Meshes from 180 to 300 elements were used 
in this analysis. Typical spacings between nodal points were 0.015 
inch radially and 0.03 inch in the z-direction. 
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Chapter III 
EXPERIMENTAL METHOD 

The objective of the experiment was to determine the extent of 
contact between two plates when bolted together. Sixteen type 304 
stainless steel plates, 4 inches in diameter, were machined to 
nominal thicknesses of 1/16, 1/8, 3/16 and 1/4 inch, 4 plates 

for each thickness. After rough machining these plates were stress 
relieved at 1875°F and ground flat to 0.0002 inch. One side of each 
plate was then lapped flat to better than one fringe of sodium light 
(11 micro-inches) in the case of the 1/8, 3/16 and 1/4 inch plates, 

and to better than two fringes in the case of the 1/16 inch plates. 
Disregarding scratches, the finish of the lapped surfaces was 5 micro- 
inches rms. Each plate had a central hole, 0.257 inch in diameter, 
for a 1/4 - 20 bolt, and two notches and two holes on the periphery 
(see Fig. 7). Two techniques were employed in determining the area of 
contact when two of these plates were bolted together. The first 
technique entailed the following procedure (see Fig. 7): 

(a) The plates were cleaned with alcohol and lens tissue. 

(b) One plate was placed on the base of the fixture shown in 
Fig. 7, lapped surface up and the two holes on the periphery 
of the plates engaged with two pins on the fixture. Spacers 
between the fixture base and plate prevented the pins from 
extending beyond the top surface of the plate. 
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(c) A second plate was placed on top of the first plate, lapped 
surfaces mating. The notches on the two plates were lined 
up with each other and with notches in the base of the 
fixture. Thus, rotation of the plates was prevented. 

(d) A standard 1/4 - 20 hex-nut with its annular bearing surface 
(0.42 inch O.D.) lapped flat was engaged on a high strength 
1/4 - 20 bolt. The nut was located about two threads away 
from the head of the bolt and served in lieu of the bolt head. 
The lapped surface of the nut faced away from the bolt head 
and since the nut was not sent home against the bolt head, the 
looseness of fit between nut and bolt offered a degree of self 
alignment. 

(e) The bolt and nut assembly described in (d) above was then 
inserted through the 1/4 inch central holes of the two plates 
and a second 1/4 - 20 lapped nut was engaged on the bolt. 

Thus the two plates were captured by the two 1/4 -20 nuts 
with the lapped surfaces of the nuts bearing against the plates. 

(f) With the torque wrench shown on the right in Fig. 7, the nuts 
were torqued down to 70 pound-inches of torque to yield a 
1100 pound force in the bolt [11]. 

(g) The position of the keys was changed to engage with only the 
lower plate and the fixture and a special spanner wrench, as 
shown in Fig. 7, was engaged with the top plate. The spanner 
wrench was restrained to move in the horizontal plane and it 
was set into motion by the screw pressing against the wrench 


handle. 
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(h) With the aid of the spanner wrench the upper plate was 

rotated relative to the lower plate several times approxi- 
mately + 5 degrees. 

Thus, the above procedure allowed for the rubbing of one plate 
relative to its mate while under a bolt force of approximately 1100 lbs. 
The remaining steps were the disassembly and the measurement of the 
extent of the contact zone which was defined by the shine due to the 
rubbing in the contact zone. It is to be noted that the boundaries of 
the contact zone as measured by the naked eye and by searching for 
marks of "polished” or "damaged" surface under a 10.5 power magnifi- 
cation are essentially the same. 

The above test was performed on 5 pairs of specimen. These were 


1. One 

0.07 

in. plate mated 

to 

a 

0.65 

in. 

plate 

2 . One 

0.126 

in. plate mated 

to 

a 

0.126 

in. 

plate 

3. One 

0.191 

in. plate mated 

to 

a 

0.192 

in. 

plate 

4 . One 

0.253 

in. plate mated 

to 

a 

0.256 

in. 

plate 

5. One 

0.124 

in. plate mated 

to 

a 

0.257 

in. 

plate 


The identical tests were repeated for 

1. One 0.124 in. plate mated to a 0.126 in. plate; and 

2. One 0.191 in. plate mated to a 0.192 in. plate, 

but in lieu of the 1/4 - 20 nuts in direct contact with the plates 
special washers, 1.000 in. O.D., 0.257 in. I.D. and 0,620 in. high, 

were interposed between the bolt head and nut. 

The diameters of the contact zones were measured with a machinist 
ruler with 100 divisions to the inch and with a Jones and Lamston 
Vertac 14 Optical Comparator. 
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The second technique used the same parts and fixture, but it 
involved autoradiographic measurements. 

Four plates, 1/4, 3/16, 1/8 and 1/16 inch thick were sent to 

Tracerlab, Inc., Waltham, Mass., for electrolytic plating with radio- 

M 

active silver Ag 110 (half life of 8 months). Each plate was 
masked except for an area on the lapped face one inch in radius. The 
plates then received a plating of copper about 5 microinches thick and 
then approximately a 5 microinch plating of silver containing the 
radioactive isotope. The resultant activity on each plate was about 
2 millicuries. 

These plates were then mated to plates of equal thickness (not 
plated) and assembled in a shielded hood as indicated in steps (a) to 

(h) above except that in the case of the pair of 1/4 inch plates care 

was taken not to rotate the plates during and after assembly and in the 
remaining cases the rotation specified in step (h) was done only once 
in one direction. 

The plates were then disassembled and the radioactive contamination 

on the plates which were in contact with the radioactive plates measured. 

The transferred activity was; 

1/4 in. plate approximately 0.05 microcuries 
3/16 in. plate approximately 3. microcuries 
1/8 in. plate approximately 0.1 microcuries 

1/16 in. plate approximately 0.4 microcuries 

It was also observed in handling that the adhesion of the silver on 
the 3/16 in. plate was poor. 
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Kodak type R single coated industrial x-ray film was then placed 
on the contaminated plates under darkroom conditions. The sensitive side 
of the film was pressed against the radioactive sides of the plates 
with a uniform load of about five pounds and left for exposure for three 
days. After three days, the film was removed and developed. The 
results are shown in Fig. 10. 
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Chapter IV 
RESULTS 

A. Pressure Distribution and Radii of Separation from Single Plate 
and Two Plate Finite Element Models. 

Using the finite element procedure described in Chapter II, the 

midplane stress distribution of single circular plates of thickness 2D, 

outer radii of 1.54 in., inner radii of 0.1 in., Poisson ratio of 0.3, 

and loaded by a constant pressure between radii A and B, Fig. 3(f), 

was computed. Computations were performed for D values of 0.1, 0.1333 

and 0.2 in. For each value of D the radius B, which defines the 

region of the symmetric external load, assumed the values of 0.31, 0.22, 

0.16 and 0.13 in. The O stress distribution at the midplane, from 

z 

the inner radius to the radius at which the above stress is no longer 
compressive, is shown in Figs. 12, 13 and 14 as a function of radius. 

The identical cases were then recomputed, using again the finite 
element method, in accordance with the two plate model shown in Figs. 4(b) 
and 5(b). These results are given in Figs. 15, 16 and 17. 

Inspection of the above figures show that the two plate model yields 
a somewhat different stress distribution in the contact zone than the 
stress distribution approximated from the single plate model, and more 
significantly, from the heat transfer point of view, tne two plate 
model yields a lower value for the radius of separation, R^, 


which 
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results in a reduction in area for heat transfer. Table 1 gives a 
comparison of the values for obtained from the two models. 

It may be observed that the single plate result of Femlund (Ref. 3, 
pp. 56, 124) is in fair agreement with the finite element results 
obtained for the single plate model. 


B. Radii of Separation from Experiment and Their Predicted Values from 
the Two Plate Finite Element Computation. 

As described in Chapter III, stainless steel circular plate 
specimen (Fig. 7) were bolted together, rotated relative to each other 
with the bolt force acting, and after disassembly the contact area of 
the joint was determined by measuring the footprints (the shiny, polished 
areas) on each plate due to the plates rubbing against each other. 
Photographs of these footprints are shown in Fig. 8. Fig. 9 also shows 
a typical footprint of the annular bearing surface of the 1/4 - 20 nut 
against a plate. All plates tested were of 304 stainless steel, 

4 inch O.D., .257 I.D., and the nominal thicknesses of the plates were 

1/16, 1/8, 3/16 and 1/4 inch. In addition to the plates fastened 
with standard nuts which gave a loading circle of radius B (Fig. 5) of 
0.211 inch, plates fastened by the special nuts described in Chpater III 
for which B was 0.5 inch V7ere also tested. 

Figure 10 shows the results of the autoradiographic tests described 
in Chapter III. For all plate pairs tested, i.e. 1/16, 1/8, 3/16 and 
1/4 inch nominal, the value of B was 0.211 inch. 

The pressure distributions and radii of separation for all the 
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above test cases were computed independently by the two plate model 
finite element analysis. Table 2 gives the test and analytical results 
for all test cases. The test results are an average of all measurements 
(minimum of six readings). A description of the analyses follows. 

Figure 18 shows the results of a two plate and a single plate model 
analysis for the 0.253 inch bolted test specimen. For Figure 19 
the external pressure distribution between radii A and B is triangu- 
lar. (The total force, however, is equal to the force exerted in the 
case of uniform pressure.) In one case, the peak external pressure is 
at A, Fig. 20(a), and in the other case at B, Fig. 20(b). Results 
of another computation which assumed a uniform displacement of 50 
microinches under each nut is shown in Fig. 21-. It is interesting to 
note that the point of separation obtained by using the two plate model 
for all variations of loading given above occurs in the range of r/A 
values of 2.73 to 2.93 while the two plate model yields separation at 
a value for r/A of 3.5. The computed deflections under the nuts are 
given in Fig. 22. 

The finite element analysis results for the 0.191 in. plate pair 
specimen are given in Fig. 23. Figures 24 and 25 show the computed 
pressure distribution and deflection patterns in the joint, respectively, 
for the 1/8 in. plate pair. In order to investigate the possible 
influence misalignments of the spanner wrench, i.e. vertical forces or 
restraints exerted at edge of plate, may have on the results of the 
experiment, the extreme case of fixing the outer edges of the plate as 
shown in Fig. 20(c) was considered. As Fig. 24 shows, within the 
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resolution of the finite element grid. size, the effect is negligible. 
This model, Fig. 20(c), and result also indicate that the influence of 
additional fasteners 2 inches away would not have an influence on the 
contact zone for the geometry considered. (However, if the distance 
between bolts is considerably reduced, then the contact area should 
increase.) The computed results for the 1/16 inch plate pair is 
given in Fig. 26. 

Figure 27 gives the finite element analysis results for the asym- 
metric case of a 1/8 in. plate bolted to a 1/4 in. plate. The model 
shown in Fig. 5(c) was used and as discussed in Chapter II, this 
model is strictly valid only if the friction in the joint prevents 
sliding between the plates. Nevertheless, the- percent discrepancy 
between the computed value and tested value (see Table 2) falls within 
the range of the symmetric cases analyzed and tested. 

In summary, the results obtained from the two plate finite element 
model and from experiment are in good agreement (Fig. 28) . 
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Chapter V 
APPLICATION 

An application of the above results for the evaluation of the 

thermal contact conductance, h , and the determination of the heat 

c 

transferred in a specific, but typical, lap joint section is illus- 
trated in this chapter. 

An aluminum lap joint in a vacuum environment, the relevant 
section and boundary conditions as shown in Fig. 29, was analyzed by 
means of a nodal analysis. The plate thickness was 0.1 in. and the 
hole diameter, 2A, was 0.2 in. The bearing surface of the bolt, 2B, 
was 0.26 in. in diameter. Because of the high conductivity and small 
thickness of the plates, no z dependence (see Fig. 29) was assumed 
for the temperature in the main body of the plate. However, heat flow 
in the z-direction in the nodes above and below the contact zone is 
considered. Qualitatively, the heat flow in the joint proceeds in the 
x-y plane from the left end (Fig. 29) toward the 0.2 in. diameter 
hole. In the vicinity of the hole, a macroscopic constriction for heat 
flow is encountered because the flow is being channeled toward the 
small contact zone. The flow of heat then encounters the microscopic 
constrictions at the contacting asperities (which determine h^) in the 
contact zone; spreads out in the x— y directions in the second plate; 
and continues to the right edge of the lap joint. 
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The material properties assumed were (refer to equation 1.1): 

H == 150,000 psi 

k = 100 Btu/hr-°F-ft (k^ = ^2 " 

a = 5.9 X 10 ^ ft. (a^ = 0^ ~ ^ ^in.) 

tan 9 = 0.1 


Assuming further, a uniform load of 46,500 psi on the loading surface 
(#10 screw; 1000 lb. bolt force) and referring to Fig. 15, curve 

g 

— = 1.3, the following interface stresses, a , contact heat transfer 
^ z 

coefficient, h , and conductance, (area)*(h ), were obtained as a 
c c 

function of inner and outer radii. (These radii define increments of 
area, the sum of which define one quarter of the contact zone.): 


outer 

r . 

inner 


h 

c 

Area x 

inch 

inch 

psi 

Btu/hr-°F-ft^ 

Btu/hr-°! 

.13 

.1 

27,900 

446,000 

16.6 

.16 

.13 

14,000 

223,000 

10.6 

.175 

.16 

3,950 

63,100 

1.7 


The conductance between nodal points were then computed and with the aid 
of the steady state heat transfer program listed in Appendix D, the 
nodal temperatures for the conditions given in Fig. 29 were computed. 

The heat transferred from the edge maintained at 20°F to the edge at 
0°F (Fig. 29) for this case was 2.88 Btu/hour. The same computation 
was repeated for the case of a bearing surface between the plate 
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and the bolt (2B) of 0.44 in. in diameter, but the bolt force was 
left unchanged. The heat transferred from the 20° F edge to the 
0°F edge in this case was 3.15 Btu/hour. In the absence of the 
joint the heat transfer along an equivalent 7 inch length of solid 
aluminum would have been 3.58 Btu/hour. This data shows that the 
thermal resistance of the contact zone (not entire 7 inch lap joint) 
was decreased from 1.52 to 0.92 °F-hr/Btu by the increase of the 
effective bolt head diameter from .26 to .44 in. It should be observed 
that the change in thermal resistance of the joint is primarily due to 
the increase in contact area and the resulting decrease in macroscopic 
constriction resistance at the hole. Also, the heat fliix in this 
example is mainly controlled by the 7 inch length and 0.1 inch 
thickness rather than the joint resistance. This emphasizes the import- 
ance of a balanced thermal design. 

For large heat fluxes where thermal strains may have an influence 
on the radii of separation, the finite element program given in 
Appendix C may be used. Also, in a non-vacuum environment the effect of 
the interstitial fluid is added in two ways. Firstly, equation (1.6) 
is applied to account for the presence of interstitial fluid in the 
contact zone, and secondly, conduction across the gaps between the plates 
and convection from the plates is considered. (Radiation heat transfer, 
if applicable, should also be included.) 
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Chapter VI 
CONCLUSIONS 

The finite element technique used in this work for the analysis 
of the pressure distribution and deformation of smooth and flat bolted 
plates under conditions of axial symmetry predicts contact areas in 
joints considerably lower than reported previously in the literature. 
These results were verified experimentally. The discrepancy between 
the previously reported results and the results reported here is due 
to the simplifying assximption made by earlier -researchers that a joint 
can be modeled as a single plate. 

The computer programs listed in the appendices will also accom- 
modate joints made up of plates of dissimilar materials and the presence 
of thermal gradients. 

Of the eleven tests performed, only one (case 3, autoradiographic) 
yielded inconsistent results. (This data point could probably be 
ignored because of the poor adhesion of the plating material which 
manifested itself by the high radioactive contamination count during 
test . ) 

The finite element analysis performed for the test specimen show 
that the gap between the 1/4 inch bolted steel specimen is 98.6 
microinches at the outer radius of the plate of 2 inches, and 1/32 
of an inch away from the radius of separation (0.35 in.), the gap is 
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only 3 microinches for the test load. This data indicates the 
difficulties previous workers have encountered in their experiments. 
(This also explains the oval shape of several of the footprints.) 
Furthermore, this data shows that the effects of surface roughness and 
the lack of flatness could have a significant effect on the size of 
contour area. 

An application of the above work to a heat transfer problem is 
illustrated in Chapter V. 
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APPENDIX A 

FINITE ELEMENT ANALYSIS OF AXISYMMETRIC SOLIDS 

The finite element method and the equations which govern the 
stresses and displacements in axisymmetric solids is given in the 
literature [9,10,12,13,15] and the procedure will be briefly summarized 
in this appendix. 

The procedure for the standard stiffness analysis method is as 
follows [15]: 

(a) The internal displacements, v, are expressed as 

{v(r,z)} = [M(r,z)] {a} (A.l) 

where M is a displacement function and a are the generalized coordi- 
nates representing the amplitudes of the displacement functions. 

(b) The nodal displacements v^ are expressed in terms of the 
generalized coordinates 

{V,} = [A] {a} <A.2) 

where A is obtained by substituting the coordinates of the nodal points 
into M. 

(c) The generalized coordinates are expressed in terms of the nodal 
displacements 

{a} = [A] ^ {v^} 


(A.3) 
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(d) The element strains, £, are evalioated 

{e} = [B(r,z)] {a} (A.4) 

where B is obtained from the appropriate differentiation of M. 

(e) The element stresses are expressed in terms of the stress- 
strain relation D 

(a(r,z)} = [D] {el = [D] [B] {a} (A. 5) 

(f) Assuming a virtual strain e and a generalized virtual 
coordinate displacement a the internal virtual work, W^, in the 
differential volumn, dV, is given by 

dW = {e}'^ (a} dV = {a}'^ [B]'^ [D] [B] {a} dV (A.6) 

1 


and the total internal virtual work is 


W. 

X 



[D] [B] 


dV ] a 


(A. 7) 


(g) The external work, W^. associated with the generalized 
displacement a is 

W = {g} (A. 8) 

e 


where ^ are generalized forces corresponding with the displacements a. 
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(h) After equating and and setting the a displacement 


to unity 


I 


{ & } = [ j [B]" [D] [B] ] a = [ k ] {a} 

Vol 


I 


where [ k ] = ) [B] [D] [B] dV 

‘‘Vol 


(A. 9) 


(A.IO) 


and which transforms to the nodal point surfaces 


k = [A ^] [ k ] [A 


(A.ll) 


(i) The stiffness matrix for the complete system is then 


n 


[K] = E [k] 

IIF=1 


m 


(A. 12) 


where n equals the number of elements and the equilibrium relationship 
becomes 

(A. 13) 


where 


{Q} 

= [K] 


n 

{Q} 

I {R} 

{r} 

= / tA-b" 


Area 


m 


(A. 14) 


(A.15) 


and P are the surface forces. 

The above procedure applies with minor modification to problems 
with thermal and body force loading. 
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The expression 

{q} = [K] {v^} (A. 16) 

represents the realtionship between all nodal point forces and all 
nodal point displacements. Mixed boundary conditions are considered 
by rewriting this equation in the partitioned form 



where v. = u. 

1 

The first part of the partitioned equation can be written as 


aa 


{u } + 

a 




and then expressed in the reduced form 


{Q*} = [K ] {u } 

aa a 

where 

{q*} = {Qj - [K^^] {u^} 


(A.19) 


(A.20) 


The matrix equation (A.19) is solved for the nodal point displacements 
by standard techniques . Once the displacement are known the strains 
are evaluated from the strain displacement relationship and the stresses 
in turn are evaluated from the stress strain relations. 

Both triangular and quadrilateral elements are used. The displace- 
ments in the r-z plane in the element are assumed to be of the form 
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+ OL^z 

(A. 21) 

+ a^r + a^z 

This linear displacement field assures continuity between elements since 
lines which are initially straight remain straight in their displaced 
position. Six equilibrium equations are developed for each triangxilar 
element . 

A quadrilateral element is composed of four triangular elements 
and ten equilibrium equations correspond to each element. 
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APPENDIX B 

finite element program for the analysis of isotropic 

ELASTIC AXISYMMETRIC PLATES (ref. 13,14) 

Input Ins t ructions : 


Sequence 

Item 



Format 

Columns 

1 

Title 



18A4 

1-72 

2 

Total number of nodal points 



15 

1-5 


Total number of elements 



15 

6-10 


Total number of materials 



15 

11-15 


Normalizing stress (NORM) 



15 

16-20 


Number of pressure cards 



15 

21-25 


(If NORM = 0, put in value of E 

in 

material 

card; 


if NORM = 1 , put in value 

^^*^vertical’ 



if NORM = -1, put in value 

^octahedral’ 



NOTE: Use NORM = 0 for this application.) 


3 

(Material property cards - one 
(a) and (b) for each material) 

set 

of 




(a) 1st card 






Material No. 



15 

1-5 


Initial a stress 

z 



FIO.O 

6-15 


Initial a stress 

r 



FIO.O 

16-25 


(b) Second Card 






E 



FIO.O 

1-10 


V 



FIO.O 

11-20 
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Card 

Sequence 


Item 


Format Column 


4 


Nodal point information (One for 
each node) 


2I5,4F10.0 


Node number 

CODE 6-10 

r-coordinate 11-20 

z-coordinate 21-30 

XR 31-40 

XZ 41-50 


If the number in column 10 is 

0 XR is the specified R-load and 
XZ is the specified Z-load 

1 XR is the specified R-displacement and 
XZ is the specified Z-load 

2 XR is the specified R-load and 

XZ is the specified Z-displacement . 

3 XR is the specified R-displacement and 
XZ is the specified Z-displacement. 


Condit ion 
free 



fixed 


Remarks 


The following restrictions are placed on the size of problems which 


can be handled by the program. 

Xtem Maximum Number 


Nodal Points 

450 

Elements 

450 

Materials 

25 

Boundary Pressure Cards 

200 




-45A- 


All loads are considered to be total forces acting on a one radian 
segment. Nodal point cards must be in numerical sequence. If cards 
are omitted, the omitted nodal points are generated at equal intervals 
along a straight line between the defined nodal points. The boundary 
code (column 10), XR and XZ are set equal to zero. 

If the number in columns 6-10 of the nodal point cards is other 
than 0, 1, 2 or 3, it is interpreted as the magnitude of an angle in 
degrees. The terms in columns 31-50 of the nodal point card are then 
interpreted as follows: 

XR is the specified load in the s-direction 

XZ is the specified displacement in the n-direction 

The angle must always be input as a negative angle and may range from 
-.001 to -180 degrees. Hence, +1.0 degree is the same as -179.0 
degrees. The displacements of these nodal points which are printed by 
the program are 

u^ = the displacement in the s-direction 
= the displacement in the n-direction 

Element cards must be in element number sequence. If element cards 
are omitted, the program automatically generates the omitted information 
by incrementing by one the preceding 1, J, K and L. The material 
identification code for the generated cards is set equal to the value 
given on the last card. The last element card must always be supplied. 

Triangular elements are also permissible; they are identified by 
repeating the last nodal point number (i.e. I, J, K, K) . 

One card for each boundary element which is subjected to a normal 
pressure is required. The boundary element must be on the left as one 
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progresses from I to J. Surface tensile force is input as a negative 
pressure. 

Printed output includes: , 

1. Reprint of input data. 

2. Nodal point displacement 

3. Stresses at the center of each element. 

Nodal point numbers must be entered cotmterclockwise around the 
element when coding element data. 

The maximum difference between the nodal point numbers on an 
element must be less than 25. However, on a nodal diagram elements 
and nodes need not be numbered sequentially. 
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n n o 


rx{N,4)=T X(N-1,4U1 
TX{N,5)=IX(N-1 

170 WRITE (6*2017) N , ( I X ( N , I ) ♦ I = 1 * 5 ) 

TF (H-N) 180*180.140 
130 IF (NUMEL-N) 300.300*130 

* s 4: « =4: ^ ^ ^ ^ ^ ^ ^ ^ 

READ AND PRINT THE PRESSURE CAROS 

300 IF(NUHPC) 290.210.290 

290 WRI TE( 6.9000 ) 

DO 200 L=1.NUHPC 

RFAD( 5.90C1 ) I BC ( L ) , J BC ( L ) . PR ( U 
200 WRI TE (6.9002 I IRC ( L) . J BC ( L ) . PR ( L ) 

210 CONTINUE 

C DETERMINE BAND WIDTH 

r^ XJ ^ 11^ ;;i; :4: ^ ;f; ijt 4: ^ ^ 3?*^ 3i« i^S ^ :j}; ^ 3^^ 5jc :?i 4s ■!?« 3Si 34 c j 5? 5j« ^ .* ^ Jj: ap 3i^: 3^ 

J = 0 

no 340 N=1,NUMEL 
on 340 1^1,4 
on 32 5 1_=1*4 
KK=IX(N,I)“IX{N*L) 

IE (KK.LT.O) KK=-KK 
TF (KK.GT^J) J=KK 
325 CONTINUE 
340 CONTINUE 

MRAN0=2^J+2 

C j^:»5jjT(c!^:^4::4e:^*35cs4c4:^«^^:4«3<«*:3r^J|«*«J!e35e3j:*:4c*W3!c**34c.3i<54{*5is*):jJt;;4J{t3jcje!^;{s*J4cj#:Ji:#**5jt3»tJjc3jc«345 4s A’:?< 

C SOLVE FOP DISPLACEMENTS AND STRESSES 

;j«; XC 4e A 4« !j: * ^ 44 ^ 3i* ^ ^ ^ 3^' ^ 3{* 3^C 5^ 43 34s 34^ 2^ ^ ^ 3?! ^ 345 * 

KSW = 0 

CALL STIFF 

IF (KSW.NE.O) 00 TO 900 

L-. 

CALL BANS CL 
/4R ITF( 6* 2052 ) 


FENT0073 

FENT0C74 

FFNT0075 

FENT0076 

FENT0077 

FENT0C78 

FENT0079 

FENT0080 

FENT0081 

FENT0G82 

FENT0083 

FENT0C84 

FENT0085 

FENT0086 

FENT0087 

FENTQ088 i 

FENT0089 

FENT0090 I 

FENT0091 

FENT0092 

FENT0093 

FENT0094 

FENT0095 

FENT0096 

FENT0097 

FENT0098 

FFNT0G99 

FFNTOIGC 

FENTOlOl 

FENT0102 

EENT0103 

FENT0104 

FENTQ1C5 

FENT0106 

FENT0IC7 

FENT0108 



NRITF (6t2025l (N,B (2^N-1),B ( 2^N ) , 1 , NUMN P ) 

C 

450 CALL STRFSSCSPLOn 

C PROCESS ALL DECKS EVEN IF ERROR 

GO TC 910 
900 WRITE (6*4C00) 

^10 WRITE 16,4001) HED 
C 

920 READ (5,1000) CHK 

IF (CHK.NE.STRS) GO TO 920 
GG TO 50 
950 CONTINUE 

WRITE (6,4002) 

CALL EXIT 

1000 FORMAT (18A4) 
loot FORMAT (1215) 

1002 FORMAT ( I5,2F10*0) 

1003 FORMAT { 2F10.0 ) 

1004 FORMAT (2F10.0) 

1005 FORMAT (3F10*0) 

1006 FORMAT (2I5,4F10*0) 

1007 FORMAT (615) 

2000 FORMAT (1HI,20A4) 

2006 FORMAT (28H0NUM8ER OF NODAL POINTS 13/ 

1 2RH NUMBER OF ELEMENTS 13) 

2007 FORMAT ( 20H0MA TE R I AL NUMBER^ 13/ 

1 25H INITIAL VERTICAL STRESS^ F10*3 ,5X, 

2 26HIN1TIAL HCRIZONTAL STRESS^ F10.3) 

2013 FORMAT (12H1N0DAL POINT ,4X, 4HTYPE ,4X, lOHR-ORD I NA T E ,4X, 
1 10H7-0RCINAT E ,10X, 6HR^LGAD ,10X, 6HZ-L0AD ) 

2014 FORMAT (I12,I8,2F14.3,2E16.5) 


FENTQ109 

FENTOllO 

FENTOlll 

FENT0112 

FFNT0113 

FFNT0114 

FENT0115 

FENT0116 

FENT0117 

FENT0118 

FENT0119 

FENT012C 

FENT0121 

FENT0122 

FENT0123 

FENT0124 I 

F6NT0125 g 

FENT0126 , 

FENT0127 

FFNTQ128 

FFNT0129 

FENT0130 

FENT0131 

FENT0132 

FENT0133 

FENT0134 

FENT0135 

FENT0136 

FENT0137 

FENT0138 

F6NT0139 

FENT0140 

FENT0141 

FENT0142 

FENT0143 

FENT0144 



P.015 FORMAT {26H0N0DAL POINT CARO ERROR N= 15) 

2016 FORMAT I49H1ELEMENT NO. I J K L MATERIAL ) 

2G17 FORMAT ( 1 1 13 ,4 16, U 12 ) 

2025 FORMAT C12H0NGDAL POINT ,6X, 1 4HR~0 I SPL AC EMENT ,6X, 14HZ-DISPLACEM 
lENT / ( I12,1P2D?0.7) ) 

2041 FORMAT (76H0M0CULUS AND YIELD STRESS NORMALIZED WITH RESPECT TO IN 
IITIAL VERTICAL STRESS ) 

2051 FORMAT! lHO,iOX ,»E» , 8X , » NU • , / , 3X , F 1 1 . 1 , FIO. 4/ ) 

2052 FORMAT! IHl) 

3G03 FORMAT U6I5) 

0 Vic ,)5< ^ 3^ ^ 3k ajc 3jf jjt # 3jj jjf sjc sj? 4 5^^ 35«3iic^Si 'ij: ;jJS5<e3^cc^ 

4000 FORMAT (//// * ABNORMAL TERM I NAT ION * > 

4001 FORMAT (//// • END OF PROBLEM * 20A4) 

4002 FORMAT (////• END OF JOB*) 

5^:4;:i^^3V:5j5 s{« J^esfc 3 ^ 5jc.3K: 3 ^ jjt 3 ^ jfejjecJcaJcjJc ;;Se3jt :^c jjc;§« 3 ^ 3jt 

9000 FnRMAT(29H0PRFSSURE BOUNDARY CONDITIONS/ 24H I J PRESSU 

IRE ) 

9001 FnRMATC2I5,FlC.O) 

9002 FnRMAT{2I6,F12.3) 

FND 

SUBROUTINE STIFF 
C 

IMPLICIT REAL*^8 (A'-H,G-Z) 

IMPLICIT INTEr3ER^2{ I-N) 

COMMON STTCP,HED(18),SIGIR(25),SIGIZ(25) , GAMMA (25 ) , ZKNOTC 25) , 

1 DEPTH{25) ,E( I 0 . 2 5 ) , $ IG ( 7 ) , R ( 450 ) , Z ( 450 ) , UR C 4 50 ) , 

2 U7 (450), ST0TAL(450,4) ,KSW 

COMMON / INTFGRZ NUMNP , NUMEL , NUMM AT , N DE PTH ,NQR M , M T YPE , I CODE (450) 
COMMON /ARG/ RRR ( 5 ) , Z ZZ ( 5 ) , S ( 1 0 , 10 ) , P ( 1 0 ) , LM ( 4 ) , DD ( 3 , 3 ) , 

1 HH( 6, 10) , RR( 4) ,ZZ(4) .0(4,4) ,H( 6, 10 ) ,Q{ 6,6) , F(6, 10) . TP ( 6) ,X 1(6) , 

2 EE( 10) , TX(450,5) 

COMMON /RANARG/ 8 ( PCO ) , A ( 900 , 54 ) , MB AND 
COMMON /PRESS/ IBC( 200 ) , J8C ( 200 ) , PR ( 200) ,NUMPC 
DIMENSION C0DE(450) 

A 3k ^ -t- * ^ ^ 4 ^ ^ ^ ^ ^ ^ ^ 4; ^ ^ ^ ^ ^ 4' )ik jjc J(t 3ik jjt ^ 
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400 CONTINUE 
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500 RETURN 
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RCC3 FORMAT ( R6H0NE GAT I V E AREA ELEMENT NO* 14) 

R004 FORMAT ( R9HOBANO WIDTH EXCEEDS ALL0WA3LE 14) 

r 4:* ^4t'!icjjS3^t^^3^^’4e«3®3$c:i!5}S3»e5jC5!c3^t;i:3!cjjc*jJt55e*#:«ti«c:S6#35t?Jt#Jit3teJ{C3i£*i«t'5««3is3jtjJt|»j>t:^3jt!jt*3itajtiJe5iC3ic:«!:#^3jc2^# 

END 

SUBROUTINE QUAD (N* VOL) 

C 

IMPLICIT REALMS (A-H,0-Z) 

IMPLICIT INTEGFR>^2< I-N ) 

COMMON STTOP, HEO( 18 ) , S I G I R ( 25 ) , S I G 1 2 < 25 ) , GAMMA { 2 5 ) t ZKNQT 125) ♦ 

1 DEPTH! 25)*EI1C,25),SIG(7 ) *R(450) , Z ( 450 ) , UR ( 450 ) ♦ 

2 UZ 1450 ) * STnTALI450*4 ) ,KSW 

COMMON /INTEGR/ NUMNP , NU MEL , NUMM AT , ND E PT H ♦NORM , MTYP E * IC00E(450) 
COMMON /ARG/ RRR ( 5 ) , ZZ Z ( 5 ) . S ( 10 ♦ 10 ) ♦ P ( 1 0 ) t LM ( 4 ) ♦ DO ( 3 , 3 ) , 

1 HH( 6*10) , RR (4 )*ZZ(4) ♦C(4,4) ♦H(6* 10) •C(6»6) ,0! 6, 10)f TP(6) f XI (6) t 

2 FEIIO) ♦ IX(450»5 ) 

COMMON /6ANARG/ B ( 900 ) ♦ A { 900 * 54 ) ♦ MBANO 

I=:IX{N,1 ) 

J = I X ( N * 2 ) 

K = T X ( N , 3 ) 

L=IX( N«4) 

C 

I 1=1 
12 = 2 
13-^3 
14 = 4 
15=5 

C DETERMINE ELASTIC CONSTANTS AMD STRESS-STRAIN RELATIONSHIP 

('_ ft /< X« . 3 C ^ 3 j{ 4 ; ■••;< * »; i' j^c 3 jc s» 3 ^ :{i ap ♦ :jc jjc ^ ;{e sje j{i jjf * ajt jJe :* :4ciat >!{ * k v >!« sfs * * * ^ ♦ * * ★ '•* sic 4c ♦ ^ ♦ :«c 

c 

CALL MPROP(N) 


FENT0239 
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FENT0309 
FENT0310 
FENT0311 
FENT0312 
FENT0313 
FFNT0314 
FENT0315 
FENT0316 
FENT0317 
FENT0318 
FENT0319 
FFNT0320 
FENT0321 
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APPENDIX C 

FINITE ELEMENT PROGRAM FOR THE ANALYSIS OF ISOTROPIC 
ELASTIC AXISYMMETRIC PLATES - THERMAL STRAINS INCLUDED (Ref. 13, 14) 


Program Capabilities : 

The following restrictions are placed on the size of problems which 
can be handled by the program. 

Item Maximum Number 


Nodal Points 

450 

Elements 

450 

Materials 

25 

Boundary Pressure Cards 

200 


Printed output includes: 

1. Reprint of Input Data 

2. Nodal Point Displacements 

3. Stresses at the center of each element. 

Input Data Format : 

A. Identification card - (18A4) 

Columns 1 to 72 of this card contain information to be printed 
with results. 

B. Control card - (5I5,F10.0) 

Columns 1-5 Number of nodal points 
6-10 Number of elements 
11 - 15 Number of different materials 



- 69 


16 - 20 Normalizing stress (see NOEM, Appendix B) 

21 - 25 Number of boundary pressure cards 

26 - 35 Reference temperature (stress free 
temperature) 

C. Material Property information 

The following group of cards must be supplied for each 
different material: 

First Card - (215, 2F10.0) 

Columns 1-5 Materials identification - any number 

from 1 to 12 . 

6-10 Number of different temperatures for 

which properties are given = 8 maximum. 

11 - 20 Initial Z stress. 

21 - 30 Initial R stress. 

Following Cards - (4F10.0) One card. for each temperature 

Columns 1-10 Temperature 

11 - 20 Modulus of elasticity - E 

21 - 30 Poisson's ratio - V 

31 - 40 Coefficient of thermal expansion 

D. Nodal Point Cards - (215, 5F10.0) 

One card for each nodal point with the following information: 

Columns 1-5 Nodal point number 

10 Number which indicates if displacements 
or forces are to be specified. 

11-20 R - ordinate 

21-30 Z - ordinate 

31-40 XR 

41 - 50- XZ 

51 - 60 Temperature 
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If the number in column 10 is 


Condition 


0 XR is the specified R-load and 
XZ is the specified Z - load 

1 XR is the specified R-displacement and 
XZ is the specified Z-load. 

2. XR is the specified R-load and 

XZ is the specified Z-displacement . 

3 XR is the specified R-displacement and 
XZ is the specified Z- displacement. 


free 



fixed 


All loads are considered to be total forces acting on a one radian 
segment. Nodal point cards must be in numerical sequence. If cards 
are omitted, the omitted nodal points are generated at equal intervals 
along a straight line between the defined nodal points. The necessary 
temperatures are determined by linear interpolation. The boundary code 
(column 10), XR and XZ are set equal to zero. 


Skew Boundaries: 

If the number in columns 5-10 of the nodal point cards is other 
than 0, 1, 2 or 3, it is interpreted as the magnitude of an angle in 
degrees. The terms in columns 31-50 of the nodal point card are then 
interpreted as follows: 

XR is the specified load in the s-direction 

XZ is the specified displacement in the n-direction 

The angle must always be input as a negative angle and may range from 

-.001 to -180 degrees. Hence, +1.0 degree is the same as -179.0 degrees. 

The displacements of these nodal points which are printed by the 


program are 
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= the displacement in the s-direction 
= the displacement in the n-direction 


E. Element Cards - (615) 


One card for each element 


Columns 1 

- 5 

Element number 

1. 

Order nodal points 

6 

- 10 

Nodal Point I 


counter-clockwise 
around element. 

11 

- 15 

Nodal Point J 



16 

- 20 

Nodal Point K 

2. 

Maximum difference 
between nodal point 

21 

- 25 

Nodal Point L 


I.D. must be less 

26 

- 30 

Material Identi- 
fication 


than 25. 


Element cards must be in element number sequence. If element cards 
are omitted, the program automatically generates the omitted information 
by incrementing by one the preceding I, J, K and L. The material iden- 
tification code for the generated cards is set equal to the value given 
on the last card. The last element card must always be supplied. 

Triangular elements are also permissible; they are identified by 
repeating the last nodal point number (i.e., I, J, K, K) . 

F. Pressure Cards - (215, IFlO.O) 

One card for each boundary element which is subjected to a 
normal pressure. 

Columns 1-5 Nodal Point I 

6-10 Nodal Point J 

11 - 20 Normal Pressure 

The boundary element must be on the left as one progresses from I to J. 
Surface tensile force is input as a negative pressure. 
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Listing ; 

Q xii2|c ;^4 jjc )}; jfe jjr. 

C FINITE ELEMENT PROGRAM FOR THE ANALYSIS OF ISOTROPIC ELASTIC 

C AXYSYMMETRIC PLATES REF FEAST 1,3 SAAS 2 

(2 * « ^ iji '4« ;^ jij .-^ jji 4; A 4; ^:feaj4 A ^ 

r. 

IMPLICIT REALMS <A”H,0-Z) 

IMPLICIT INTEGER^2( I-N) 

COMMON STTOPtHEDI 18),SIGIR(25),SIGIZ(25) , GAMMA (25) ,ZKNOT (25) , 

1 OEPTH (25 ) , E( 8*4*25 ) t S TG( 7 )* R(450) * Z ( 450 > ,UR ( 4 50 ) , T T ( 3) ♦ 

2 U7(450)*STOTAL(450*4) * 

3 T( 450 ) .TFMP*Q*KSW 

COMMON /INTFGR/ NUMNP, NUMEL , NUMMAT* NDEPTH,NORM, MTYPE , ICODE(450) 
COMMON /ARC/ R RR ( 5 ) , Z Z Z ( 5 ) , S ( 10 , 10 ) , P ( 1 0 ) » LM( 4 ) * DO ( 3 , 3 ) , 

1 HH (6 ,10 ) ,RR(4) ,ZZ(4) ,C(4*4) ,H(6, 10) *D(6*6) , F(6»10 ) tTP( 6) tXI (6) , 

2 EE{10),IX{450,5) 

COMMON /BANARG/ R ( 900 ) , A ( 900 , 54 ) , MB AN D 
CGMMGN/PRESS/ IBC(200 ) , J0C ( 200 ) , PR ( 200),NUMPC 
DATA STRS /» ak "sjc 3f< 3j« I f 

A is ,'k 5k ^ 4: 3?t aja jjt 3fe 3^ ^ * j£ 3jc jjc j^e ^ )^c s^c 3^ :^s 3j< ^ 

C READ AND PRINT CONTROL INFORMATION 

^ *it # Jjc 3{C )Jc ap 

50 READ { 5,1000,ENO=950) HED 
WRITE (6,2000) HEO 
C 

READ( 5, 1001 ) NUMNP, NU MEL , NUM MAT . NORM , NUMPC , 0 
WR J TE(6, 2C06 ) NUMNP , NOME L ,NUMMAT , NUMPC , 0 
IF (NORM) 65,65,66 
66 WRITE (6,2041) 

READ AND PRINT MATERIAL PROPERTIES 

65 CONTINUE 
C 

DO eC M=1,NUMMAT 

READ (5, 1012) MTYPE, NUMTC, SIG I Z ( MTYPE ) ,SIGIR( MTYPE) 

WRITE (6, 2011 )MTYPF,NUMTC, SIG IZ( MTYPE) ,STGIR( MTYPE) 


EEWTOCCl 

FEWT0002 

FEWT0C03 

FEWT0004 

FFWT0005 

FEWT0006 

FEWT0007 

FEWT0008 

FEWTOC09 

EEWTOO 10 

FEWTOOIl 

EEWT0012 

FEWT0013 

FEWT0014 

FEWT0015 

FEWT0016 I 

FEWT0017 -j 

FEWT0018 

FEWT0C19 ' 

FFWT002C 

FEWT0021 

FEWT0022 

FEWT0023 

FEWT0024 

FEWT0025 

FEWT0026 

FEWT0027 

FEWT0C28 

FEWT0029 

FFWT0030 

FEWT0C3I 

FEWT0032 

FEWT0033 

FEWT0034 

FEWT0035 

FEWTQ036 
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R FAD (5, 10 11 ) ( (E( I, J^i^TYPF), J=l,4), I = ltNUMTC) FEWT0G37 
WRITE (6,2010) ((E( I,J,MTYPF) , J=1 , 4 ) , 1= 1 , NUMTC ) FEWT0C38 
DO 81 I=NUMTC,8 FEWT0039 
00 81 J=l,4 FEWT0040 
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130 READ (5,1CC7) M , ( I X { M , I) , I =1 , 5 ) 

140 N=N+1 

IF (M-N) 170, 170 ,150 
150 IX(N, I )=IX(N-1,1) + 1 
TX<N,2) = IX( N-1 ,2 ) +1 
TX{N,3)=IX(N-1,3)+1 
IX(N,4)^ IX(N~1 ,41+1 
IXCN,5)=IXCN~1 ,5) 

170 WRITE (6,?C17) N , ( I X { N , I) , I = 1 , 5 ) 

IF IM-NI 180,180,140 
180 IF (NUFEL-N) 300,300,130 

C READ AND PRINT THE PRESSURE CAROS 

300 IF(NUHPC) 290,210,290 
290 WR ITEC 6,9000 ) 

DO 200 L^l.NUMPC 

READ! 5,9001 ) I 8C ( L ) , J BC (L ) , PR( L I 
200 WRITE (6,9002 ) IBC( LI , JBC ( L) , PR ( L ) 

210 CONTINUE 

* S5S * ^ 3^j j§t :$t 3^ 3^ ^ ^ Jjc 3^ * Jjc ijs 34: 3{{ jjt ^ }jc 3^c ^ 1!^ ^ * Jjc j|c ^ ;{c 2^ 35s ^ 5jc ft 

DETERMINE BAND WIDTH 
J=0 

no 340 N=1,NUHEL 
00 340 1=1,4 
DO 325 1=1,4 
KK=IX(N, I )~TX( N,U 
IF (KK.LT.OI KK=-KK 
IE (KK.GT.J) J=KK 
325 CONTINUE 
340 CONTINUE 

MBANn = 2’!^ J+2 

C SOLVE FOR DISPLACEMENTS AND STRESSES 


FEWT0073 

FEWT0074 

FEWT0075 

FEWT0076 

FEWT0077 

FEWT0C78 

FEWT0C79 

FEWT0080 

FEWT0081 

FEWT0082 

FEWT0C83 

FEWT0084 

FEWT0C85 

FEWT0086 

FEWT0C87 

FEWT0088 

FEWT0089 

FEWT0090 

FEWT0091 

FEWT0092 

FEWT0093 

FEWT0094 

F6WT0C95 

FEWT0096 

FEWT0C97 

FEWT0098 

FEWTCC99 

FEWTOIQO 

FEWTOlOl 

EEWT01C2 

FEWT01C3 

FEWT0104 

FEWT0105 

FEWT0106 

FEWT0I07 

FEWT01C8 
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KSW = 0 

CALL STIFF 

IF (KSW.NE^O) GO TO 900 
C 

CALL BANSOL 
yRrTF(6,2052 ) 

WRITE ( 6*2025) (N*B (2^N-U,B ( 2^N ) , 1 * P ) 

C 

450 CALL STRESS! SPLOT) 

C PROCESS ALL DECKS EVEN IF FRROR 

GO TO 910 

900 WRITF (6*4000) 

910 WRT TF { 6,4001 ) HED 

920 READ (5,1000) CHK 

IF (CHK. NE, SIRS) GO TO 920 
GO TO 50 
950 CONTINUE 

WRITE (6,4002) 

CALL EXIT 

;^t *1}; :4rjj< A ^ w ;^e j): ^ 3|' ^ sje x« ^ 3 ^ j(c ^ jJ; 3 ^ ^ Asjt; jje 3^s jjtaft jjc 

)Ji; i^s A Jiii i; j{t ^ ^ jjs V ajc 4c S[c 

ICOO FORMAT (18A4) 

1001 F0RMAT(5I 5,F10.0) 

1 002 ^0R^4AT ( I5,2F10.0) 

1003 FnRMAT(2F10.0) 

1004 FORMAT (2F10.0) 

1005 FORMAT (3F10.0) 

1006 FORMAT! 21 5,5F1C,0) 

1007 FORMAT (615) 

1011 FORMAT (4F1 0.0) 

1012 FORMAT! 21 5,2F1C.C) 

4:4^4£4cA^4c4£*4:3^«4:4e4l4t4t*;4^4tx«4!4:^^t4cJ{t4c4c4e4s5^t4cJjcJ{^4c4f4tl^^5;c*4;4c4cJ^c4c4t4:’?«^5C^3^c3^c3i^^4‘5^'*4c4:‘§fs^«♦4^’^^^ 

2000 FORMAT (1H1,20A4) 


FEWT01C9 
FEWTOllO 
FEWTOlll 
FEWT0112 
FEWT0113 
FEWT0114 
FEWT0115 
FEWT0116 
FEWT0117 
FEWT0118 
FEWT0119 
FEWT012C 
FEWT0121 
FEWT0122 
FEWT0123 
FEWT0124 
FEWT0125 
FEWT0126 
FEWT0127 
FEWT0128 
FEWT0129 
FEWT0130 
FEWTG131 
FEWT0132 
FEWT0133 
FEWT0134 
FEWT0135 
FEWT0136 
FEWT0137 
FEWT0138 
FEWT0139 
FEW! 01 40 
FEWT0141 
FEWT0142 
FFWT0143 
FFWT0144 



?006 FORMAT (//, 

1 30H0 NUMBER OF NODAL POINTS 13 / 

2 30H0 NUMBER OF ELEMENTS 13 / 

3 30H0 NUMBER OF DIFF. MATERIALS 13 / 

4 30H0 NUMBER OF PRESSURE CARDS 13 / 

5 30H0 REFERENCE TEMPERATURE E12*4) 

2010 FORMAT (15H0 TEMPERATURE 15X SHE 15X 6HNU 15X 6HALPHA 9X 

1/4F20.8) 

2011 FORMAT 1 1 THOM ATER 1 AL NUMRER= 13* 30H, NUMBER OF TEMPERATURE CAROS = 

1 13,25H INITIAL VERTICAL STRESS^ F10.3*5X, 

2 27H INITIAL HORIZONTAL STRESS^ F10*3I 

2013 format I12H1N00AL POINT ,4X, 4HTYPE *4X* lOHR-ORO I NA TE *4X* 

1 10H7~nRCINATE * 1 OX , 6HR-L0AD *10X* 6HZ-L0 AD * lOX , 4HTEMP ) 

2014 FORMAT 1112, I 8*2F14. 3, 2E16.5*F14* 3) 

2015 FORMAT (26H0N0CAL POINT CARD ERROR N= IS) 

2016 FORMAT (49H1ELEMENT NO. I J K L MATERIAL ) 

2017 FORMAT (1113,416,1112) 

2025 FORMAT (12H0N0DAL POINT *6X, 14HR-DI SPLACEMENT *6X, 14HZ-0ISPLACEM 
lENT / ( I 12, IP2D20.7) ) 

2041 FORMAT {76H0M0CULUS AND YIELD STRESS NORMALIZED WITH RESPECT TO IN 
IITIAL VERTICAL STRESS ) 

2051 FOPMAT( 1HO,10X*»E» ,8X, *NU’ ,/ ,3X,F11. 1,F10.4/ ) 

2052 FORMAT(lHl) 

3003 FORMAT (1615) 

^ XS i; * ^ T«c # 3«e ^ jiK aje ^ Jjt 5? * 

4000 format (//// * ABNORMAL TERMINATION*) 

4001 FORMAT (//// » END OF PROBLEM • 20A4) 

4002 FORMAT {////* END OF JGB«) 

^ ^ ^ * 3^ 3^ 1(1 3jc sjt :4c sk Jlc S5C 3^ 5^c ^ >is « sit - t s?s « ^ ^ 3^ >se * jje * * 3^t * t 

9000 FORMAT (29H0PRESSURE BCUNDARY CONDITIONS/ 24H I J PRESSU 

IRE ) 

9001 FGRMAT(2I5,F10.0) 

9002 F0RMAT(2I6,F12.3 ) 

FND 

SUBROUTINE STIFF 


FEWT0I45 

FEWT0146 

FEWT0147 

FEWT0148 

FEWT0149 

FEWT015Q 

FEWT0151 

FEWT0152 

FEWT0153 

FEWT0154 

FEWT0155 

FEWT0156 

FEWT0157 

FEWTOlSa 

FEWT0159 

FEWT0160 

FEWT0161 

FEWT0162 

FEWT0163 

FEWT0164 

FEWT0165 

FEWT0166 

FEWT0167 

FEWT0168 

FEWT0169 

FEWT0170 

FEWT0171 

FEWT0172 

FEWT0173 

FEWT0174 

FEWT0175 

FEWT0176 

FEWT0177 

FEWT0178 

FEWT0179 

FEWT0180 



■") n o 


C 


IMPLICIT REALMS <A~H,C-7) 

IMPLICIT INTeGER^2( r-N) 

COMMON STTOP*HED< 18)tSIGIR(25)»SIGIZ(25) ,GAMMA(25) *ZKNCT (25) , 

1DEPTHC2 5) ,E{ 8^4,25 } ♦ S IG ( 7 ) , R ( 450 ) » Z ( 450 ) » UR ( 4 50 ) t TT( 3 ) , 

2 UZ(450),STOTAL(450*4) , 

3 T(450) ,TEMP*0,KSW 

COMMON /INTEGR/ NUMNP, NUMEL , NUMMAT * NOEPTH,NORM, MTYPE , I CODE (450) 
COMMON /ARG/ RRR ( 5 ) .ZZZ ( 5 ) , S ( 10» 10 ) » P ( 10) ,LM( 4 > »00< 3 ,3 ) » 

1 HH(6,10) t RR(4)tZZ(4) ,C<4,4) ,H( 6,10) ,0(6,6) ,F(6, 10) f TP(6) ,XI (6) , 

2 EE(IO) ,IX(450,5) 

COMMON /BANARG/ B ( 900 ) , A ( 900 ,54 ) , M8AN0 
COMMON /PR ESS/ IBC ( 200 ) , JBC ( 200 ) , PR ( 200 ) , NUMPC 
DIMENSION C0DE(450) 

-;3e Jjt Jjc 3^c ^ :sis ^ ^ 3^ 35t' jfe * *!c Jje sSc ^ ^ 3^ 3^ ^ sjt ^ 

INI TI AL IZATION 

NB=^27 
N0 = 2=^NB 
N02=2=S'NUMNP 
DO 50 N=1,N02 
B(N)=0.0 
DO 50 M=1,ND 
50 A(N,M)=0.0 

C FORM STIFFNESS MATRIX 

on 210 N=1,NUMFL 
C 
C 

90 CALL OUAD(N,VOL) 

IF (VOL) 142,142,144 
14? WRITE (6,2003) N 
KSW=1 
GO TO 210 
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CnSA=:OCOS( CODEC I ) ) 

29? BUJ-1) = B{ JJ-l H?.X=^(C0SA^DZ + SINA^DRJ 
R(JJ)=BCJJHZX’?'(S INAX^DZ-COSA^OR ) 

3C0 CONTINUE 
310 CONTINUE 
C DISPLACEMENT B.C. 

C 

00 400 H=1*NUMNP 
U=UR(M) 

KX=TC0DE(M)+1 

CO TO (400,370,390,380 ,KX 
370 CALL M0DIFYCN,U,ND2) 

GO TO 400 

330 CALL MQDIFYC N,U,ND2) 

390 U=UZ(M) 

N = N + 1 

CALL MODIFY(N,U,N02) 

400 CONTINUE 
C 

500 RETURN 

2003 FORMAT ( 26H0NEGAT I VE AREA ELEMENT NO. 14) 

2004 FORMAT (29H0RAND WIDTH EXCEEDS ALLOWABLE 14) 

END 

SUBROUTINE OUAC(N,VGL) 

C 

IMPLICIT REAL^^S (A-H,0-Z) 

IMPLICIT INTEGER=^2( I-N) 

COMMON STTOP,HED( 18) , SI G I R ( 2 5 ) , S I G I Z ( 25 ) , GAMMA (25) , ZKNOT (25) , 

1DFPTH(25),E(8,4,25),STG(7),R(450) ,Z(45C) , UR (4 50) ,TT(3) , 

2 UZ(450) ,STOTAL(450,4) , 

3 T( 450 ), TEMP, Q,KSW 

COMMON / INTEGR/ NUMNP, NUMEL , NUMMAT , N DEPTH, NORM, HTYPE , I CODE ( 450 ) 
COMMON /ARG/ R R R ( 5 ) , Z Z Z ( 5 ) , S ( 1 0 , 10 ) , P ( 1 0 ) , LM ( 4 ) ♦ DO ( 3 , 3 ) , 


FEWT0289 

FEWT0290 

EEWT0291 

FEWT0292 

FEWT0293 

FEWT0294 

FEWT0295 

FEWT0296 

FEWT0297 

FEWT0298 

FEWT0299 

EEWT03C0 

FEWT0301 

FEWT0302 

FEWT0303 

FEWT0304 ' 

FEWT03C5 o 

EeWT0306 , 

FEWT0307 

EEWT0308 

FEWT0309 

FEWT0310 

FEWT0311 

FEWT0312 

FEWT0313 

FEWT0314 

FEWT0315 

EEWT0316 

FEWT0317 

FEWT0318 

FEWT0319 

FEWT0320 

FEWT0321 

FEWT0322 

FEWT0323 

FEWT0324 



n o 


C 


c 


c 


c 

c 

c 


c 

c 


1 HH(6,10) *RR( 4 ),7.Z(4) »C(4t4) *H( 6 t 10) f Di 6t6) tF{ 6f IG) f TP(6) »XI (6) f 

? EEdO) , IX(450,5) 

/BANARG/ B ( 900 ) t A ( 900 ♦ 54 ) , HB AND 

I = !X(N,1 ) 

J= I X ( N • 2 ) 

K=IX{N,3) 

L=IX(N,4) 


103 

104 


TO 

71 

105 


11 = 1 
T2 = 2 
13 = 3 
14=4 
15^5 

THERMAL STRESSES 
TEMP={ T( n+T( J )+T(K)+T(L) )/4.0 


DO 103 M=2.8 

TF( E( M, 1 . MTYPE )-TEMP) 103,104^104 


CONTINUE 


RATI0=0.0 

nFN=E( M, 1 ,MTYPE)-E(M-1 ,1 ,MTYPE) 
IF(DEN) 70,71,70 

RAT 10= ( TFMP-e ( M-1 , 1 , MTYPE ) )/DEN 


no 

FF 


105 KK=1,3 ^ 

( KK ) = E( M-1 ,KK+1, M TYPE )+R ATI 0^< E ( M , KK + 1 , MTYPE ) ~E ( M- 1 , KK + 1 ♦ MT YPE 


) ) 


TFMP=TEMP-Q 

nETERMINE ELASTIC CONSTANTS AND ST RESS-STR A I N RELATIONSHIP 


CALL MPKCP(N) 


88 DO 110 M=l,3 

110 TT(M)-(C{M,1)<-C(M,2)+C(M,3))=^^EE(3)^^TEMP 
FORM GUADRILATERAL STIFFNESS MATRIX 
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APPENDIX D 

STEADY STATE HEAT TRANSFER PROGRAM FOR BOLTED JOINT 


Program Capacity i 50 nodal points 

Output Data: 

(a) Input data 

(b) Inverse of matrix 

(c) Nodal temperature 

(d) Given and calculated augmenting vector 
and residual error 

Input Data Sequence ; 

A. Case identification (12A4) followed by two blank cards 

B. Card (II) with a 1 

C. Card (17) with dimension of matrix 

D. Card (II) with a 1 

E. Cards (II, 3(213, E15.8)) with node indices started in the 
first 13 field followed by conductance between these nodes. 
Only input from lower node number to higher node number 
required (since the conductance from node i to j equals 
the conductance from j to i.) Each card has three groups 
of z node numbers followed by a conductance value except 
the last card. Last card could have 1, 2 or 3 groups and 
has a 1 in column 1. 

F. Cards (II, 3(l6, E15.8) with number of node followed by 
conductance from the node to ground node which is at 
specified temperature. Each card has 3 groups of node 
number followed by conductance. The II field is skipped 
except for the last card for ground conductances which can 
have 1, 2 or 3 fields and the first column has a 1. A 
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node can be connected to only one ground node. 

G. Same as F above, but code temperature specified for ground 
node instead of the conductance value. 

H. Same as F above, but code internal power dissipation for the 
particular node instead of the conductance value. 



STFAHY STATE HEAT TRANSFER PROGRAM BOLTED JOINT SSHTOOGl 

DIME^SICN IQENT (12 ), A 1050,050 ),AA(050, 050 ) ,8{ 50),8I{ 50 )t SSHT0C02 

1DC( 50)tRES( 50),ACON( 50),TACCN( 50),Q( 50) SSHT0003 

WRITE (6, 23) SSHT0004 
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(see Figs. 12 - 17) 






















TABLE 2 


Test and Analytical Results for Radii of Separation of Bolted Plates (see Fig. 5) 




2B 

in. 

Separation Diameters, 2 in. 

% Discrepancy Between 
Computed Values and 
Tested Values 

Case 

■L/ 

in. 

"Rubbing 

Test" 

Autoradiographic Test 

Computed 




Range 

Average 

Range 

Average 


Rub. Test 

Autorad. Test 

1 

.065 

,422 

.42-. 48 

.45 

. 41-. 46 

.44 

.488 

7.8 

9.8 

2 

.124 

.422 

.50-. 53 

.51 

.4 - .6 

.55 

.554 

7.9 

.7 

3 

.191 

.422 

.58-. 64 

.62 

.76-. 81 

* 

.78 

.620 

0 

25.8 

4 

,253 

.422 

0 

1 

.72 

.68-. 73 

irk 

.7 

.700 

2.9 

0 

5. 

Unmat cl 
ed Pail 
.124/ 
.257 

1- 

^ .422 

.54-. 58 

.56 

— 

— 

.588 

4.8 

— 

6. 

.124 

1.0 

1.06-1.10 

1.09 

— 

— 

1.104 

1.3 

— 

7. 

.191 

1.0 

1.11-1.17 

1.16 

— 

— 

1.210 

4.1 

— 


Original x-ray film shows hole in plate and 0.6 inch diameter zone more distinctly than remainder of 
area sensitized by the radioactive contamination. Loose radiographic contamination observed during 
test , 

Assembled and disassembled radioactive and non-radioactive plates without rotating plates relative to 
each other. 
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FIG. 2. ROETSCHER's RULE OF THUMB FOR PRESSURE 
DISTRIBUTION IN A BOLTED JOINT 
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s- 



104 



FIG. 3(b) 
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(a) Actual Plates 


(b) Finite Element 
Ideal izati on 



(c) Single Annular Ring 
Element 



FIG. 4. FINITE ELEMENT IDEALIZATION OF TWO PLATES IN CONTACT 



z 


(b) Finite Element Model for Plates of Equal Thickness 


z 



(c) Finite Element Model for Plates of Unequal Thickness 


FIG. 5. FINITE ELEMENT MODELS 
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(a) Plates Intersect, R too small 

o 


I 



FIG. 6, EXAMPLES OF UNACCEPTABLE SOLUTIONS 




FIG. 7. PLATE SPECIMEN, BOLT AND NUTS, FIXTURE AND TOOLS. 
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FIG. 8(a). FOOTPRINTS ON MATED PAIR OF 1/16 INCH PLATES. 
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FIG. 8(b). FOOTPRINTS ON MATED PAIR OF 1/8 INCH PLATES. 
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FIG. 8(c). FOOTPRINTS ON MATED PAIR OF 3/16 INCH PLATES. 





FIG. 8(d). FOOTPRINTS ON MATED PAIR OF 1/4 INCH PLATES. 



5 - 



FIG. 8(e). FOOTPRINTS ON MATED PAIR OF 1/8 AND 1/4 INCH PLATES. 

FIG. 8. FOOTPRINTS ON THE MATING SURFACES OF 1/16 - 1/16, 

1/8 - 1/8, 3/16 - 3/16, 1/4 - 1/4, and 1/8 - 1/4 
PAIRS. (A = .128, B = .21) 
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FIG. 9. FOOTPRINT OF NUT ON PLATE. 
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FIG. 10 (a). 1/16 INCH 

PAIR 



FIG. 10(b). 1/8 INCH 

PAIR 
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FIG. 10(d). 


1/4 INCH 
PAIR 



FIG. 10. X-RAY PHOTOGRAPHS OF CONTAMINATION TRANSFERRED 

FROM RADIOACTIVE PLATE TO MATED PLATE. 1/16, 1/4, 
3/16, 1/4 INCH PAIRS. (A = .128 1n., B = .21 in.) 
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FIG. n. FREE BODY DIAGRAM FOR TWO PLATES IN CONTACT. 





















FIG. 16. INTERFACE PRESSURE DISTRIBUTION IN A BOLTED JOINT 
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r/A 


FIG. 18. 


FINITE ELEMENT ANALYSIS RESULTS FOR 1/4 INCH PLATE PAIR. 
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FIG. 19. PRESSURE IN JOINT, TRIANGULAR LOADING 
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FIG. 21. PRESSURE IN JOINT, UNIFORM DISPLACEMENT UNDER NUT. 



DEFLECTION - INCHES 














FIG. 23. FINITE ELEMENT ANALYSIS RESULTS FOR 3/16 INCH 

PLATE PAIR. 








FIG. 24. FINITE ELEMENT ANALYSIS RESULTS FOR 1/8 INCH PLATE PAIR. 






DEFECTION OF ONE PLATE AT JOINT 


L33 



— FINITE 


FIG. 25. GAP DEFORMATION FOR FREE AND FIXED EDGES 
ELEMENT ANALYSIS, 1/8 INCH PLATE PAIR. 








FIG. 26 
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FIG. 27. FINITE ELEMENT ANALYSIS RESULTS FOR 1/8 INCH PLATE 
MATED WITH 1/4 INCH PLATE. 
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FIG. 28. COMPARISON BETWEEN TESTED AND MEASURED SEPARATION RADII. 






